Electrostatic interactions between graphene layers and their environment. 
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We analyze the electrostatic interactions between a single graphene layer and a Si02 susbtrate, 
and other materials which may exist in its environment. We obtain that the leading effects arise 
from the polar modes at the SiC>2 surface, and water molecules, which may form layers between the 
graphene sheet and the substrate. The strength of the interactions implies that graphene is pinned 
to the substrate at distances greater than a few lattice spacings. The implications for graphene 
nanoelectromechanical systems, and for the interaction between graphene and a STM tip are also 
considered. 
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I. INTRODUCTION. 

Graphene is a versatile two dimensional material whose 
singular electronic and mechanical properties show a 
great potential for applications in nanoelectronicspii 2 -' 3 - 
Since free floating graphene is subject to crumpling^ the 
presence of a substrate, and the environment that comes 
with it, is fundamental for its stabilization. Hence, this 
environment will have direct impact in the physical prop- 
erties of graphene. Though the influence of the substrate 
and other elements of the surroundings has been taken 
into account in different ways in the literature, the exact 
part that these are playing is not yet fully understood. 

On the one hand, the differences observed between 
samples grown on different substrates constitute an open 
issue. Most experiments have been carried out in 
graphene samples deposited over Si02, or grown over SiC 
substrates^ and a better understanding of how graphene 
properties are expected to change would be worthy. On 
the other hand, there is the question of characterizing all 
the effects that a particular environment has on graphene 
electronic and structural properties. 

Concerning electronic properties, it has been suggested 
that the low temperature mobility of the carriers is deter- 
mined by scattering with charged impurities in the Si02 
substrate^ and the effect of these charges can be sig- 
nificantly modified by the presence of water molecules. 8 
Actually, the very polar modes of Si02 give a good 
description of the finite temperature corrections to the 
mobility ^ 10 ' 11 Supporting this idea, recent experiments 
show that graphene suspended above the substrate has a 
higher mobility i 12 ' 13 

Experiments also seem to reveal a very important 
role played by the substrate in the structural proper- 
ties of graphene. STM measurements suggest that sin- 
gle layer graphene follows the corrugations of the Si02 
substrate ; 14 ' 15 and experiments on graphene nanoelec- 
tromechanical systems (NEMS) indicate that the sub- 
strate induces significant stresses in few layer graphene 
samples.™ Moreover, the interaction between graphene 
and the substrate determines the frequency of the out of 
plane (flexural) vibrations, which can influence the trans- 
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FIG. 1: a) Sketch of the system studied in the text. Inter- 
action effects: b) Interaction with water molecules attached 
to hydroxyl radicals at the substrate, c) Interaction with po- 
lar modes at the surface of the substrate, d) van der Waals 
interaction between the graphene sheet and the metallic gate. 



port properties at finite temperatures i 17 ' 18 

In order to shed light on the influence of the environ- 
ment on the graphene properties, we analyze the char- 
acteristic energies of interaction with the substrate and 
other materials present in the experimental setup. This 
allows us to evaluate the relative importance of the differ- 
ent interactions in the binding and mechanical response 
of the graphene layer. We also provide estimates of quan- 
tities such as equilibrium distances, typical lengthscales 
of corrugations, and frequencies of vibration, which can 
be measured in principle in current experimental setups. 

Throughout the paper we concentrate on Si02, though 
results are easily generalized to other substrates. Partic- 
ularly, we consider: i) the van der Waals forces between 
graphene and the metallic gate below the Si02 substrate, 
ii) the electrostatic forces between the graphene layer 
and the polar modes of the substrate, iii) the electro- 
static forces between graphene and charged impurities 
which may be present within the substrate and iv) the 
electrostatic forces between graphene and a water layer 
which may lay between graphene and the substrate^ A 
sketch of the setup studied, and the different interaction 
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mechanisms, is shown in Fig[TJ We will also mention the 
possibility of weak chemical bonds between the graphene 
layer and molecules adjacent to it , 19 ' 20 although they will 
not be analyzed in detail. We do not consider a possible 
chemical modification of the graphene layer; 21 ' 22 which 
would change its transport properties. 

The general features of the electrostatic interactions to 
be studied are discussed in the next section. Then, we an- 
alyze, case by case, the different interactions between the 
graphene layer and the materials in its environment. Sec- 
tion IV discusses the main implications for the structure 
and dynamics of graphene, with applications to graphene 
NEMS and the interaction between graphene and a STM 
tip. The last section presents the main highlights of our 
work. 



II. ELECTROSTATIC INTERACTIONS 
BETWEEN A GRAPHENE LAYER AND ITS 
ENVIRONMENT. 

The electrons in the tt and tt* bands of graphene are po- 
larized by electromagnetic potentials arising from charges 
surrounding it. The van der Waals interactions between 
metallic systems, and metals and graphene can be ex- 
pressed as integrals over the dynamic polarizability of 
both systems. Those, in turn, can be written in terms 
of the zero point energy of the plasmons i 23 ' 24 The in- 
teraction between the graphene layer and a polarizable 
dielectric like SiC>2 is also given by an integral of the po- 
larizability of the graphene layer times the polarizability 
of the dielectric. The latter can be approximated by the 
propagator of the polar modes, which play a similar role 
to the plasmons in a metal. The interaction between the 
graphene and static charges of electric dipoles depends 
only on the static polarizability^ 

We will calculate these interactions using second order 
perturbation theory, assuming a perfect graphene sheet 
so that the momentum parallel to it is conserved. The 
corresponding diagrams are given in Fig. [2j All interac- 
tions depend, to this order, linearly on the polarizabil- 
ity of the graphene layer. In ordinary metallic systems, 
the Coulomb interaction is changed qualitatively when 
screening by the graphene electrons is taken into account 
through an RPA summation of diagrams. This is not the 
case for undoped graphene. There, the Random Phase 
Approximation leads to a finite correction ire 2 /8hv-p ~ 1 
to the dielectric constant, which does not change signifi- 
cantly the estimates obtained using second order pertur- 
bation theory. 

The response function of a graphene layer at half filling 



is 
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where N s = N v — 2 are the valley and spin degeneracy. 
This expression is obtained assuming a linear dispersion 
around the K and K' points of the Brillouin Zone. It 



is valid up to a cutoff in momentum A ~ a -1 and en- 
ergy to c ~ v-pA, where a is the lattice spacing. Beyond 
this scale, the susceptibility has a more complex form, 
and it is influenced by the trigonal warping of the bands. 
The component of the electrostatic potential induced by 
a system at distance z from the graphene layer with 
momentum q is suppressed by a factor e~^ z . Hence, 
the integrations over q can be restricted to the region 
< q = \q\ ^ q-max ~ z~ x . The combination of a term 
proportional to e _ '^ z and scale invariant quantities such 
as the susceptibility in Eq. |T]) leads to interaction en- 
ergies which depend as a power law on z. In general, 
we will consider only the leading term, neglecting higher 
order corrections^ 

The calculation described above, which is valid for a 
single graphene layer at half filling, can be extended to 
other fillings and to systems with more than one layer. 
In all cases, the calculations are formally the same, and 
the interaction energies can be written as integrals over 
energies and momenta of the susceptibility of the sys- 
tem being considered, which replaces the susceptibility 
of a single layer, Eq. |T]). The susceptibility of a doped 
single layer is well approximated by that of an undoped 
system, Eq. (flj, for momenta such that q > Anal- 
ogously, the susceptibilities of a stack of decoupled layers 
of graphene and multilayered graphene become similar 
for q > t±/hvp2^, where t± is the hopping in the per- 
pendicular direction. The susceptibility of a single un- 
doped plane of graphene, Eq.([T]) is an increasing func- 
tion of q, so that the integrals are dominated by the 
region q — q max ~ z~ x . Hence, if q max > fc F or 
Qmax 3> t±/hvF, the interaction energies do not change 
appreciably from the estimates obtained for a single layer. 
The corrections can be obtained as an expansion in pow- 
ers of either k F z, or (t±z)/hvp. The expression Eq. (TT|) 
can therefore be considered as the lowest order expan- 
sion in these parameters. For z ~ lnm, t± ~ 0.35eV 
and carrier densities such that n ~ 10 10 — 10 12 cm~ 2 , we 
obtain k F z ~ 10" 2 - 10" 1 and t ± /hv F ~ 1CT 2 - 10" 1 . 
In the following, we will analyze mostly the interaction 
energies using the expression in Eq. (TT]) for the graphene 
polarizability. 



III. INTERACTIONS WITH SPECIFIC 
ENVIRONMENTS. 



A. Metallic gate. 

We describe the metallic gate as doped Si, separated 
from the graphene layer by a 300nm thick slab of SiC>2 
dielectric. For the Si doping and voltages applied, most 
of the charge in the Si gate is concentrated on a layer 
of about 10 nm thickness^ much narrower than the dis- 
tance to the graphene sheet, so that the gate is effectively 
two dimensional. We describe the susceptibility of the 



3 






x 



FIG. 2: Lowest order diagrams which contribute to the in- 
teraction between: a) graphene and a metal, b) graphene 
and a polar dielectric, and c) graphene and a static charge 
distribution. The thin red bubble stands for the graphene 
susceptibility. The thick green bubble represents the metallic 
susceptibility. The wavy green line stands for the propagator 
of a phonon mode in the dielectric. Crosses stand for static 
charge distributions, and dashed lines represent the electro- 
static potential. 



gate as that of a dirty two dimensional electron gas: 

dn Dq 2 



Xgate(q,iu) = - 



dfi Dq 2 



(2) 



where D = VF gate l ga te is the diffusion coefficient of the 
electrons in the gate, VF gate is the Fermi velocity, l ga te is 
the mean free path, and dn/ da is the bare compressibil- 
ity, given by the density of states at the Fermi level (see, 
for instance, Ref. l3ll ). 

The interaction between the graphene layer and the 
gate is given by: 



V q( Z ) 



(3) 



being e the static dielectric constant of the SiC"2 sub- 
strate. 

The lowest order contribution to the energy in pertur- 
bation theory has the following form: 
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gate 
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7^V 2 (z)xG{q,iu)Xgate(q,iu). (4) 



For future reference, note that we use the symbol E for 
energies per unit area, and £ for total (integrated) ener- 
gies. The resulting integrals can be calculated analyti- 
cally in the limit z s = D /Avp <C z: 
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The dependence on z 3 log(z/2; s ) was obtained in Ref. 

a 

We take, as representative parameters for the gate and 
the graphene layer, D s» 10 -3 m 2 /s, vf = 10 6 m/s, z = 
300 nm, dn/du ~ g(E F ) = 0.04 eV" 1 A" 2 and e = 4 for 
the SiC>2 substrate. These parameters lead to interaction 
energies of order ~ 10~ 8 meV A -2 . 



B. Polar dielectric. 



electric fields induced by the surface polar modes of 



SiO- 



32,33,34.35 



The coupling can be written as: 



(6) 



H I = Y J M qPq (&, + &!,), 
i 

where p q is the electron density operator and b q ,b q the 
creation/destruction operators for phonons, and M 2 = 
(h 2 Vp)ge~ 2qz I \qa) is the interaction matrix element, 
with g a dimensionless coupling constant. In Si02 we 
have two dominant phonon modes at KQi = 59 meV and 
htt 2 = 155 meV, with gi = 5.4 • 10~ 3 and g 2 = 3.5 • 10~ 2 
respectivelyi 1 ^ 

The lowest order contribution to the energy is given 
by: 

^1=EE / ^XG{qM\M q {z)?D?\qM, (7) 



where we have introduced the free phonon propagators: 

(8) 



D\ '(q,iuj) = - 
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The calculation can be again carried out analytically. 
In the limit z <C h = Wp/f2j we obtain: 
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which gives a z~ 2 dependence on the distance. In the op- 
posite limit z I j, which can be of interest in suspended 
graphene experiments, one obtains: 



p(2) = 
subs 
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Let us give some numerical estimates for both expres- 
sions. In the case of graphene deposited over the sub- 
strate, we have z ~ 1 nm (see Ref. [2]) and U -C z, having 



interaction energies of order E 



(2) 



-4xlO" 1 meV A" 2 . 



For suspended graphene, z ~ 300 nm, and the energies 
are of order E^ bs ~ — 10 _8 meV A -2 , i.e., of the same 
order than the contribution from the gate. 



C. Charges within the substrate. 

In this case the calculations are done considering that 
effectively all the charge is concentrated close to the sur- 
face of the SiC>2 dielectric. The second order correction 
to the energy, averaged over the charge distribution, is: 



(11) 



The interaction between the graphene layer and 
the SiC>2 substrate can be expressed in terms of the 



where we consider a Coulomb interaction v q between 
graphene electrons and charges that is statically screened 
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by the effective dielectric constant at the interface, (e + 
l)/2. Again, this contribution can be carried out analyt- 
ically: 

2 1 

Ttn-imp 1 (12) 



P 2 - 



2e 2 



1 



2?iwp 2z 



This interaction has a z~ 4 dependence, as the image po- 
tential in ordinary metals. In this case, however, this 
behavior arises from the combination of a vanishing den- 
sity of states and lack of screening in graphcnc. 

Reasonable values for the impurity concentration in 
graphene are in the range ni mp ~ 10 10 — 10 12 cm" 2 ^ 
Setting z ~ lnm, typical interaction energies are of the 

order E ch 10~ 4 - 1(T 2 meV A~ 2 . 

The present result is only valid for graphene samples 
close to the substrate. If the distance to the latter is 
larger than the typical distance between charges, d imp ~ 
^/riimp ~ 1 — 10nm, the electrons feel the net effect of the 
effective charge in the substrate. This is zero in average, 
as there should be a compensated number of positive and 
negative charges. However, if we consider a finite region 
of the substrate, fluctuations can locally give rise to a 
net effective charge. This can be estimated by replacing 
imp = n lmp l 2 — > y/ n imp l 2 in the total energy £ ch 



V 



E^l 2 , where I is the lateral sample size. 



D. Layer of water molecules. 

The properties of the SiC>2 surface are dominated, for 
thermally grown SiC>2 layers, by the presence of abun- 
dant silanol (SiOH) groups^ whose surface density is 
about 5 x 10 14 cm -2 , unless extra steps like thermal an- 
nealing in high vacuum are taken during the fabrication 
process! 37 ! 38 ' 39 Silanol sites are active centers for water 
absorption, so that the SiC>2 surface becomes hydrated 
under normal conditions, 40 which is probably the case of 
most of the graphene samples produced by mechanical 
cleavage 4 Moreover, several layers of water may cover 
the SiC>2 surface, lying between the oxide surface and 
the graphene samples after the graphene deposition. An 
analogous situation has been shown to happen in exper- 
iments with carbon nanotubes deposited on SiC>2- 41 

The water molecule has an electric dipole, p w = 6.2 x 
1CP 30 C m « 0.04e nm. Typical fields applied in present 
experimental setups are £ ~ 0.1 V ran" 1 . The energy 
of a water dipole when it is aligned with this field is 
4 meV ~ 50K, so that, at low temperatures, it will be 
oriented along the field direction, perpendicular to the 
substrate and the graphene layer. For this reason, in the 
following we assume that the water molecules are not 
charged, and their dipoles are aligned perpendicular to 
the substrate and the graphene layer. This arrangement 
can be considered an upper bound to the interaction en- 
ergy with a neutral water layer, as inhomogeneities and 
thermal fluctuations will induce deviations in the orien- 
tation of the dipoles, and will lower the interaction en- 
ergy. Note, however, that for high applied electric fields, 



a charging of water molecules of the order Qh 2 o ~ 0.1 |e| 
has been reportedi 41 - The presence of these extra charges 
would considerably enhance the interaction between the 
graphcnc layer and the water molecules. 

A water molecule which is located at a distance z from 
the graphene layer induces an electrostatic potential: 



z) = 2tt Pv 



(13) 



This potential polarizes the graphene layer and gives rise 
to an interaction energy in a similar way to the static 
charges discussed in the preceding section. The lowest 
order contribution to the energy is: 



E 



(2) 



(ep w ) 
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(14) 



where n w is the concentration of water molecules and the 
z~ 3 behavior arises from the dipolar nature of the inter- 
actions. For z — 0.3nm, which is the approximate thick- 
ness of a water monolayer ; 42 ' 43 the interaction energy is 

(2) 

^ water ~ — I2n w meV which, for a typical water concen- 
tration n w = 10 15 cm~ 2 , yields E^1 ter ~ 1 meV/A 2 . 



The expression in Eq. ([14]) can be extended to a semi- 
infinite stack of water layers. For simplicity we take a 
distance z between graphene and the uppermost layer of 
water molecules equal to the interlayer distance. In this 
case we obtain: 
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< ^2 77 n 



C(3) 



6 hv F {2zf 



(15) 



where £(3) rs 1.202 is Riemann's zeta function. The 
present result indicates that the first water layer is the 
one that mostly contributes to the binding. 



E. Van der Waals interaction between graphene 
layers. 

For comparison, in this Section we evaluate the van 
der Waals interaction between two graphene layers at the 
equilibrium distance. Using the same approximations as 
for the other contributions, we recover the result of Ref. 

m 



E, 



(2) 

G-G 
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16hvF (2z 



»3' 



(16) 



For z = 0.3nm, this expression gives an interaction en- 
ergy of 30 meV A~ 2 . This estimate is similar to other 
experimental and theoretical values of the graphenc- 
graphene interaction ) 44 ' 45 and is at least one order of 
magnitude greater than the other contributions analyzed 
earlier. 
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IV. ANALYSIS OF THE RESULTS. 

A. Comparison of the different interactions. 

Numerical estimates for the different interaction en- 
ergies obtained for reasonable values of the parameters 
are listed in Table HI The present results show that the 
leading interactions are those between graphene and the 
polar modes of the SiC>2 substrate, and between graphene 
and a possible water layer on top of the substrate. Both 
effects are of similar order of magnitude in the present 
approximation, where we have assumed that the water 
molecules are aligned in the direction normal to the sub- 
strate. 

The interactions for multi-layer graphene samples can 
be obtained by adding the separate contributions from 
each layer. The different dependences on distance imply 
that the relative strength of the interactions in samples 
with many layers can change compared to the results of 
Table I. For instance, the effects of the polar substrate 
oc z~ 2 and of charged impurities oc z^ 1 , which are of 
longer range, sum up more effectively than the binding 
effect of water: the z~ 3 decay of the graphene-water in- 
teraction suggests that only the first graphene layer is 
affected by the presence of water on the substrate. For 
the same reason, the presence of several layers of aligned 
water molecules should not increase the binding, as only 
the closest layer effectively contributes to the interaction 
energy. On the other hand, the binding effect of water 
could be enhanced if the molecules were allowed to ro- 
tate freely, therefore approaching the high polarizability 
of liquid water jS or if they were partly ionized by the ap- 
plied field^i leading to additional charges similar to the 
Coulomb impurities present in the SiC>2 substrate. 

It should be noted that we have considered here only 
long-range electrostatic interactions, for which reliable 
expressions can be obtained, in terms of well understood 
material parameters, like the molecular polarizability, 
electric dipoles, or surface modes. Still, there is a signif- 
icant uncertainty in some parameters, like the distance 
of the relevant charges to the graphene layer and the 
concentration of charged impurities and water molecules. 
We have not analyzed the possible formation of chemi- 
cal bonds between the carbon atoms and the water or 
silanol groups at the SiC-2 surface. Calculations based 
on the Local Density Functional approximatio n 19 ' 20 sug- 
gest that individual molecules can (weakly) bind to a 
graphene layer with energies of 10 — 50 meV, although 
it is unclear how these estimates are changed when the 
molecules interact at the same time with the graphene 
layer and the substrate. 



B. Corrugation of the graphene layer induced by 
the substrate. 



The attractive forces calculated in the preceding sec- 
tion imply that graphene is bound to the Si02 substrate, 
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TABLE I: Interaction energy per unit area for the mecha- 
nisms studied in the paper. For the numerical estimates we 



have used typical concentrations of 10 10 — 
impurities and 10 15 cm -2 water molecules 



10 cm charged 



as observed in experiments. Our previous analysis does 
not include the short range repulsive forces which de- 
termine the equilibrium distance. We assume that the 
total energy near the surface is the sum of the terms an- 
alyzed above, which have a power law dependence on the 
distance, and a repulsive term, E rep {z) = e rep {z^ / z n ), 
which also decays as a power law at long distances, with 
zq an undetermined length scale. For simplicity, we 
assume that the leading attractive term is due to the 
presence of a water layer, which behaves as E water = 
—€ w (zq/z 3 ). The total energy per unit area is thus: 



E(z) = 



-rep 



At the equilibrium distance, z eqi we have: 
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so that: 
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(19) 



Hence, the order of magnitude of the pinning potential 
induced by the environment on the out of plane modes of 
graphene is given by AT oc E water (z e q) / z^ q ~ 10~ 2 - 10~ 4 
meV A -4 . Defining the out of plane displacement as 
h(r), the energy stored in a corrugated graphene layer is: 



£ w 1 / d 2 r [n(Ah) 2 + Kh 2 ] , 



(20) 



where k « leV is the bending rigidity of graphene . 46 ' 47 
For modulations h(r) defined by a length scale I, the 
bending energy dominates if I <C I* = (k/AT) 1 / 4 , while 
the graphene layer can be considered rigidly pinned to 
the substrate if I 3> I*. Using our previous estimates, we 
find I* ~ lOA, so that the graphene layer should follow 
closely the corrugations of the substrate. 
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FIG. 3: Sketch of the deformation of a nanoelectromechanical 
device studied in the text. 



The pinning by the substrate implies that the disper- 
sion of the flexural modes becomes: 




(21) 



where p is the mass density of the graphene layer. At 
long wavelengths, lim/j^o Wk = <^o ~ 10~ 4 — 10 _3 meV ~ 
- 10~ 2 K. 

The estimates obtained above also allow us to analyze 
the bending of graphene NEMS due to the interaction 
with the material below, at distance <ii 16 ' 48 We assume 
that the lateral dimension of the graphene cantilever is 
I, and the maximum displacement of the graphene layer 
from a flat position is h. A sketch of the graphene can- 
tilever is shown in Figj3] We consider the force induced 
by charged impurities in the substrate below the can- 
tilever, as this is the contribution which decays more 
slowly with the distance to the graphene layer (cf. Table 
HI . If the distance of the cantilever to the substrate is d, 
and supposing d^> h, the gain in energy due to the defor- 
mation of the graphene layer is AS ~ e c hlzoh/ y/rii mp d 2 . 
We have again defined z and e c h by rescaling E c h{d) = 
e ch z o/d, with z nu 1 nm and e c /, ~ 10~ 4 — 10~ 2 meV 
A~ 2 depending on the density of impurities. The factor 



is included, as already mentioned at the end of 



Section III.C, to take into account the effect of having 
an overall neutral distribution of charge in the substrate. 
This energy should compensate the elastic response to 
the deformation, A£ e i ~ nh 2 /l 2 , leading to an equilib- 
rium value: 
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(22) 



Reminding the dependence of e c h on the density of impu- 
rities, Eq. (|22|) results in an overall behavior h cx y/ni mp . 
For structures such that d ~ 300 nm, one finds that the 
suspended graphene sheet deforms significantly (i.e. h 
becomes comparable with d) for lenghts greater than a 
few /im. 

In the case of very pure substrates (i.e. neglecting 
the presence of charged impurities), the attractive force 
that bends the graphene sheet would be determined 
by the next corrections to the energy, which decay as 
z~ 3 . Rewriting i? (2) — (e w + e ph + e G )(z /z) 3 with 
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FIG. 4: Estimate of the threshold voltage as function of 
graphene-tip separation needed to detach a graphene layer 
from the substrate. The inset shows a sketch of the geometry 
considered in the text. 



e p h — cq ~ O.lmeV/A 2 and e w ~ O.OlmeV/ A 2 , we see 
that those interactions are dominated by the coupling to 
the gate and to the polar modes, which are of compa- 
rable magnitude. A calculation similar to the one pre- 
sented above would yield deformations of the order of 
h ~ 1A x (l/d) A 1 resulting even in these cases to unstable 
graphene sheets for lenghts greater than a few fim. 



C. Interaction with a metallic tip in an STM 
experiment 

It is known that STM tips on graphite surfaces some- 
times deform the surface graphene layer ; 49 ' 50 and the un- 
derstanding of these deformations can be of interest for 
current research on graphene i 14 ' 15 ' 51 The analysis of the 
electrostatic interactions between a graphene layer and 
its environment allows us to estimate possible deforma- 
tions induced by an STM tip. We analyze the setup 
sketched in the inset of Fig01 The tip has lateral dimen- 
sion I and it is located at a distance d from a graphene 
layer. This graphene layer interacts with an underly- 
ing substrate, and a voltage V is applied between the 
graphene layer and the tip. We consider three interac- 
tions: 

i) An attraction between the tip and the graphene 
layer, which tends to deform the graphene, in the way 
shown in Fig[4j We assume that this energy is purely 
electrostatic. A simple estimate can be obtained by de- 
scribing the setup as a capacitor where the area of the 
plates is I 2 , the distance between the plates is d, and the 
applied voltage is V. The interaction energy is of the 
order: 



yap 

tG-U P « 8ne2d , 
where we define V in energy units. 



(23) 
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ii) The pinning of the graphene layer to the substrate. 
This contribution opposes the deformation of the layer. 
We write it as: 



pin 



e ■ I 2 



(24) 



where e p i n is the pinning energy per unit area. As typi- 
cal values, we will use 1 meV/A 2 for graphene on a water 
layer, and 30 meV/A 2 for graphene interacting with an- 
other graphene layer, as in graphite. 

iii) The rigidity of the layer against flexural deforma- 
tions. This term tends to keep the layer flat. The de- 
formed region is likely to be (at least) as large as the 
size of the STM tip. As a result, an upper bound to the 
energy stored in a deformation is: 



I 2 



The graphene layer will be deformed when: 



&G—tip ^ £pin ~\~ £el- 



(25) 



(26) 



Note that the approximations involved in obtaining the 
various terms are valid only if d > a. 



We consider a situation where e 



pin j 1 



and I are fixed. 
Eq. ([26|) implies that the layer is deformed if the voltage 
exceeds a threshold: 



V>V th (d) W J87T 



ne 2 d 3 



(27) 



can be obtained in terms of a few well understood micro- 
scopic parameters, and they have a simple dependence 
on the distance between the graphene layer and the sys- 
tem which induces the electrostatic field. The analysis 
presented here should give reliable estimates of the or- 
der of magnitude of the different binding energies, and 
of their relative strength. We have not considered the 
possible formation of chemical bonds, which may alter 
the results when the distances between the carbon atoms 
in the graphene layer and the surrounding materials is 
sufficiently small. 

We find that the leading effects arise from the polar 
modes of the Si02 substrate, and water which may form 
layers on top of it. A summary of the main results is pre- 
sented in Table [U The interaction energies with systems 
with N layers can be obtained, to a first approximation, 
by adding the contributions from each layer. The esti- 
mated magnitude of the interactions suggests that a sin- 
gle graphene layer is pinned to the substrate on length 
scales greater than a few lattice spacings, ~ lOA. The 
electrostatic binding modifies the long wavelength, out 
of plane flexural modes, which acquire a finite frequency, 
cj ~ 10~ 4 — 10~ 3 meV. The long range forces consid- 
ered here can also induce large deformations in graphene 
NEMS. Besides, we have analyzed the possibility of de- 
formations of the graphene layer by an STM tip. We find 
that a voltage drop of 0.5 - 2 V between the tip and the 
sample at distances 1 - 10 A is sufficient to deform the 
graphene layer. 



Assuming I ~ 10a and d ~ a, we see that the domi- 
nant contribution comes from the pinning term Eq. (|24p . 
Hence, in the physically relevant range a < d <C I, we can 
write: 



V th (d) 



8ire v 



(28) 



independent on the tip size. The threshold values for 
graphene on SiC-2 are of about 0.5-2 V for d ~ 1 - 10A, 
as schematically shown in Fig[4l 



V. CONCLUSIONS. 

We have analyzed the electrostatic interactions be- 
tween a graphene layer and the polarizable materials 
which may be present in its environment, for samples 
deposited on SiC>2. The strength of these interactions 
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